disp('tjalff');

node_name = 'n2';
n = 2;

Num_point = 400;
mumu = 1.0;

H = 0.01;
vT =0.001:0.001:0.04;
vP = zeros(1, length(vT));

for iT = 1:length(vT)
    
    T = vT(iT);
  
    [vk, vep, ck, cX, cu] = fn_solve_tba(T,H,mumu,n,Num_point);
    [vk_dH, vep_dH, ck_dH, cu_dH] ...
        = fn_solve_tba_dH(T, H, vk, vep, ck, cu);

    pressure_dH = fn_int_pressure_dH(vk, vep, vep_dH, T);
    vP(iT) = pressure_dH;
    
    fprintf('finished T = %f \n', T);
    fprintf('==========================\n\n');
    fprintf('==========================\n\n');
    
end

eptc = {vT,vP};

mkdir(['/home/yicong/matlab/bosons/datas/',node_name]);
save(['/home/yicong/matlab/bosons/datas/',node_name,'/eptc.mat']','eptc','-v7.3');

disp('tjsjlftj');